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Abstract 

Stochastic dynamical systems arise as models for fluid particle motion in geophysical flows with 
random velocity fields. Escape probability (from a fluid domain) and mean residence time (in a fluid 
domain) quantify fluid transport between flow regimes of different characteristic motion. 

We consider a quasigeostrophic meandering jet model with random perturbations. This jet is 
parameterized by the parameter (3 = 2^ cos(#), where fi is the rotation rate of the earth, r the earth's 
radius and 9 the latitude. Note that Q. and r are fixed, so f3 is a monotonic decreasing function of the 
latitude. The unperturbed jet (for < /3 < |) consists of a basic flow with attached eddies. With 
random perturbations, there is fluid exchange between regimes of different characteristic motion. We 
quantify the exchange by escape probability and mean residence time. 

For an eddy, the average escape probability for fluid particles (initially inside the eddy) escape into 
the exterior retrograde region is smaller than escape into the jet core for < (3 < 0.3333, while for 
0.3333 < P < |, the opposite holds. 

For a unit jet core near the jet troughs, the average escape probability for fluid particles (initially 
inside the jet core) escape into the northern recirculating region is greater than escape into the southern 
recirculating region for < (3 < 0.115, while for 0.385 < (3 < |, the opposite holds. Moreover, for 
0.115 < (3 < 0.385, fluid particles are about equally likely to escape into either recirculating regions. 
Furthermore, for a unit jet core near the jet crests, the situation is the opposite as for near the jet 
troughs. 

The maximal mean residence time of fluid particles initially in an eddy increases as (3 increases from 
to 0.432 (or as latitude decreases accordingly), then decreases as (3 increases from 0.432 to § (or as 
latitude decreases accordingly). However, the maximal mean residence time of fluid particles initially 
in a unit jet core always increases as (3 increases (or as latitude decreases). 
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1 Stochastic dynamics: Escape probability and mean residence 
time 



Stochastic dynamical systems are used as models for various scientific and engineering problems. We 
consider the following class of stochastic dynamical systems 

x = ai{x,y) +bi(x,y)wi, (f) 
V = a 2 {x,y) + b 2 (x,y)w 2 , (2) 

where wi(t), 102 (i) are two real independent Wiener processes (white noises) and 01,02,61,62 are given 
deterministic functions. More complicated stochastic systems also occur in applications ([J], [25], \W\). 

EH. m, G3) 

For a planar bounded domain D, we can consider the exit problem of random solution trajectories of 
Q-© from D. To this end, let dD denote the boundary of D and let T be a part of the boundary dD. 
The escape probability p(x, y) is the probability that the trajectory of a particle starting at (x, y) in D first 
hits dD (or escapes from D) at some point in T, and p{x,y) is known to satisfy (JH|> and references 
therein) 



-y i (x,y)p xx + -bl(x,y)pyy + a 1 (x,y)p x + a 2 {x,y)py = 0, (3) 



\b\{x,y)p xx + h 

fir = 1, (4) 
p\od-t = 0. (5) 



Suppose that initial conditions (or initial particles) are uniformly distributed over D. The average escape 
probability P that a trajectory will leave D along the subboundary T, before leaving the rest of the 
boundary, is given by (e.g., E2) 

P = J J^p(x,y)dxdy, (6) 

where |D| is the area of domain D. 

The residence time of a particle initially at (x, y) inside D is the time until the particle first hits dD (or 
escapes from D). The mean residence time u(x,y) is given by (e.g., |32, [21], HH| an d references therein) 



-bl(x,y)u x 



-b 2 (x,y)u yy 



ai(x,y)u x 



a 2 (x,y)u y 
u\dD 



(7) 

(8) 



2 A quasigeostrophic jet model 

The Lagrangian view of fluid motion is particularly important in geophysical flows since only Lagrangian 
data can be obtained in many situations. It is essential to understand fluid particle trajectories in many 
fluid problems. Escape probability (from a fluid domain) and mean residence time (in a fluid domain) 
quantify fluid transport between flow regimes of different characteristic motion. Deterministic quantities 
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like escape probability and mean residence time can be computed by solving Fokkcr-Planck type partial 
differential equations. 

We now use the these ideas in the investigation of meandering oceanic jets. Meandering oceanic jets 
such as the Gulf Stream are strong currents dividing different bodies of water. 

Recently, del-castillo-Negrete and Morrison ([H]), and Pratt et al. [22]) have studied models for 

oceanic jets. These models are dynamically consistent to within a linear approximation, i.e., the potential 
vorticity is approximately conserved. Del-castillo-Negrete and Morrison's model consists of the basic flow 
plus time-periodic linear neutral modes. 

In this paper, we consider an oceanic jet consisting of the basic flow as in del-castillo-Negrete and 
Morrison ([E]), plus random-in-time noise. This model incorporates small-scale oceanic motions such as 
the molecular diffusion (d>|), which is an important factor in the Gulf Stream (UJ, d|). The irregularity 
of RAFOS floats (|3J, [33], [201) also suggests the inclusion of random effects in Gulf Stream modeling. 

This random jet may also be viewed as satisfying, approximately in the spirit of del-castillo-Negrete 
and Morrison (0), the randomly wind forced quasigeostrophic model. Several authors have considered the 
randomly wind forced quasigeostrophic model in order to incorporate the impact of uncertain geophysical 
forces ([30], H3], \Ri\, EHl! [H; 0])- They studied statistical issues such as estimating correlation coeffi- 
cients for the linearized quasigeostrophic equation with random forcing. There is also recent work which 
investigates the impact of the uncertainty of the ocean bottom topography on quasigeostrophic dynamics 

(CHI). 

The randomly forced quasigeostrophic equation takes the form ([221) 

A^ + J(^,A^) + ^ iC = ^, (9) 

where W(x,y, t) is a space-time Wiener process (white noise). The stream function would have a random 
or noise component (C2j, [HI)- Note that (5 is the meridional derivative of the Coriolis parameter [2Ej El- 
i.e., = 2j2 cos(0), where £1 is the rotation rate of the earth, r the earth's radius and 9 the latitude. Since 
Q and r are fixed, (3 is a monotonic decreasing function of the latitude. 

The deterministic meandering jet derived in del-castillo-Negrete and Morrison (0) is 

^(x, y) = — tanh(y) + a sech 2 (y) cos(kx) + cy, 

where 

a = 0.01, c = ~(1 + J/J), k = ^2(lT^ = p), < (3 < |. 

This ty(x,y) is an approximate solution of the usual quasigeostrophic model 

A^t + J(i>, A^j) + /3rl) x - . (10) 

With random wind forcing or molecular diffusive forcing in the stochastic quasigeostrophic model ©, 
the stream function would have a random or noise component. We approximate this noise component by 
adding a noise term to the above deterministic stream function ty(x,y), that is, in the rest of this paper, 
we consider the following random stream function as a model for a quasigeostrophic meandering jet, 

^(x, y) = — tanh(y) + a sech 2 (y) cos(kx) + cy + noise. 
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Figure 1: Unperturbed jet: e = and (3 = 1/3. 

The equations of motion for fluid particles in this jet then have noise terms. We further approximate 
them as white noises (or Wiener processes) 

dx = -^ y dt + y/edwi, (11) 

dy = ^> x dt + ^fedw 2 , (12) 

or more specifically, 

dx = [sech 2 (y) + 2a sech 2 (y) tanh(y) cos(fcir) — c]dt + \/edwi, (13) 

dy = — ak sech 2 (y) sin(kx)dt + \fedw2, (14) 

where < e < 1, and w\ (t) , wi (t) are two real independent Wiener processes (in time only). The calcula- 
tions below are for e = 0.001. Note that j3 is now the only parameter 1131 ITU as a is given and c, and k 
depend only on (3 € [0, |]. 

When e = 0, the deterministic jet consists of the jet core and two rows of recirculating eddies, which are 
called the northern and southern recirculating regions. Outside the recirculating regions are the exterior 
retrograde regions; see Figure 

3 Escape probability 

We take D to be either an eddy or a piece of jet core (see Figures El EJ- This piece of jet core has the 
same horizontal length scale as an eddy, and it is one period of the deterministic jet core (note that the 
deterministic velocity field is periodic in x). Thus we call this piece of jet core a unit jet core. 
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Figure 2: An eddy: /3 = 1/3 
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Figure 3: A unit jet core near a trough: (3 = 1/3 
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Figure 4: Escape probability of fluid particles (initially in an eddy in Figure |5J) exiting into the exterior 
retrograde region: /3 = 1/3 

From 0, (0J, I©, the escape probability p(x,y) that a fluid particle, initially at {x,y), crosses the 
subboundary T of the domain D satisfies 

eAp + [sech 2 (y) + 2a sech 2 (y) tanh(j/) cos(fca:) — c]p x — ak sech 2 (y) sin(kx)p y — 0, (15) 

P\r = 1, (16) 
p\ dD -r = 0. (17) 

We take T to be either top or bottom boundary of an eddy or a unit jet core (see Figures GJ- We 
numerically solve this elliptic system for various values of (3 between and |. In the unit jet core case, we 
take periodic boundary condition in horizontal (meridional) x direction, with period =£■. 

A piecewise linear, finite element approximation scheme was used for the numerical solutions of the 
escape probability p(x, y), and the mean residence time u(x, y), described by the elliptic equations and 
(0 , respectively. Using a collection of points lying on the boundary, piecewise cubic splines were constructed 
to define the boundary of the eddy, and the top and bottom boundaries of the jet core. Computational 
(triangular) grids for the eddy and the jet core were then obtained by deforming regular grids constructed 
for an ellipse and a rectangular region, respectively. The computed escape probability crossing the upper 
or lower boundary of an eddy or a unit jet core, for the case of f3 = |, are shown in Figures 0] El an d 

Suppose that the fluid particles are initially uniformly distributed in D (an eddy or a unit jet core). We 
can also compute the average escape probability P that a particle will leave D along the upper or lower 
subboundary T, using the formula ©; see Figures [H] and |5J 

For an eddy, the average escape probability for fluid particles (initially inside the eddy) escape into 
the exterior retrograde region is smaller than escape into the jet core for < (3 < 0.3333, while for 
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Figure 5: Escape probability of fluid particles (initially in an eddy in Figure EJ) exiting into the jet core: 
P = 1/3 




Figure 6: Escape probability of fluid particles (initially in the unit jet core in Figure exiting into the 
northern recirculating region: = 1/3 
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Figure 7: Escape probability of fluid particles (initially in the unit jet core in Figure 0) exiting into the 
southern recirculating region: (3 = 1/3 

0.3333 < j3 < |, the opposite holds (Figure |SJ). Thus (3 = 0.3333 is a bifurcation point. Also, the average 
escape probability for fluid particles escape into the exterior retrograde region increases as (3 increases from 
to 0.54 (or as latitude decreases accordingly), and then decreases as (3 increases from 0.54 to | (or as 
latitude decreases accordingly). Thus (3 = 0.54 is another bifurcation point. The opposite holds for the 
average escape probability for fluid particles escape into the jet core. 

For a unit jet core near the jet troughs, the average escape probability for fluid particles (initially 
inside the jet core) escape into the northern recirculating region is greater than escape into the southern 
recirculating region for < (3 < 0.115, while for 0.385 < (3 < |, the opposite holds. Moreover, for 
0.115 < (3 < 0.385, fluid particles are about equally likely to escape into either recirculating regions 
(Figure EJl. 

Furthermore, for a unit jet core near the jet crests, the situation is the opposite as for near the jet 
troughs. 

4 Mean residence time 

The mean residence time u(x, y) of a fluid particle, initially at (x, y) in either an eddy or a piece of jet core 
(see Figures 12 13J), satisfies 

eAu + [sech 2 (y) + 2a sech 2 (y) tanh(y) cos(fcx) — c]u x — ak sech 2 (y) sm(kx)u y — —1, (18) 

u\ dD = (19) 

The mean residence times of fluid particles in an eddy or a unit jet core are shown, for the case of 
(3 = \, in Figures HOI U 
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Figure 8: Average escape probability for a fluid particle in an eddy: — for exiting into the exterior 
retrograde region; for exiting into the jet core region. 
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Figure 9: Average escape probability for a fluid particle in the unit jet core near a trough: — for exiting 
into the northern recirculating region; for exiting into the southern recirculating region. 
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Figure 10: Mean residence time in an eddy: (3 
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Figure 11: Mean residence time in a unit jet core: 
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Figure 12: Maximal value of mean residence time, as a function of (3, in an eddy. 

The maximal mean residence time of fluid particles initially in an eddy increases as (3 increases from 
to 0.432 (or as latitude decreases accordingly), then decreases as (3 increases from 0.432 to | (or as latitude 
decreases accordingly); see Figure IT2*1 However, the maximal mean residence time of fluid particles initially 
in a unit jet core always increases as (3 increases (or as latitude decreases accordingly); see Figure ITSI 

5 Discussions 

The present work on fluid particle motion in random flows takes into account of fluid particle diffusive as 
well as advective motion. There has been recent work on fluid particle advective motion (molecular diffusion 
ignored) in time-periodic, quasi-periodic and aperiodic flows (periodic — > quasi-periodic — > aperiodic — * 
random); see, for example, |2*3|. [33]. |3l"j. [2]. [£]. [2T] and |2*2*]. 

Our work on random particle motion does not require that the random part is small. However it 
does require that the flow can be decomposed into steady or unsteady deterministic (drift) and random 
(diffusion) parts; Otherwise the Fokker-Planck formalism does not hold. 
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Figure 13: Maximal value of mean residence time, as a function of (3, in a unit jet core. 
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